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Abstract 

The equilibrium and stability of the coupled elastic 
lead/lag, flap, and torsion motion of a cantilever rotor blade 
in hover are addressed, and the influence of several higher- 
order terms in the equations of motion of the blade is deter- 
mined for a range of values of collective pitch. The blade 
is assumed to be untwisted and to have uniform proper- 
ties along its span. In addition, chordwise offsets between 
its elastic, tension, mass, and aerodynamic centers are as- 
sumed to be negligible for simplicity. The aerodynamic 
forces acting on the blade are modeled using a quasi-steady, 
strip-theory approximation. 


1. Introduction 

An important problem in helicopter dynamics is the 
determination of the dynamic response and aeroelastic sta- 
bility associated with the rotor blades. Considerable atten- 
tion has been directed to rotary-wing aeroelasticity prob- 
lems, and it is now widely recognized that such problems 
are inherently nonlinear. Hodges and Dowell 1 developed 
a comprehensive set of differential equations of motion, 
with quadratic nonlinearities, describing the flap-lead/lag- 
torsional dynamics of slender, rotating extensional rotor 
blades undergoing moderately large elastic deformations. 
An ordering scheme based on a small parameter e was in- 
troduced in Ref. 1 to systematically neglect higher-order 
terms in the equations. Some important linear terms of 
order e 3 were kept in the equations such as aerodynamic 
damping terms in the lead/lag and torsional differential 
equations and inertia terms in the torsional differential 
equation. Nonlinear terms of 0(e 3 ) were systematically ne- 
glected. The equations of motion developed in Ref. 1 were 
used in Ref. 2 to investigate the stability of the elastic mo- 
tion of a uniform cantilever rotor blade in the hover flight 
condition. 

A set of 0(e 3 ) nonlinear differential equations describ- 
ing the flexural-flexural-torsional motion of inextensional 
beams undergoing moderately large deformations was de- 
rived by Crespo da Silva and Glynn and used by the same 
authors to analyze the response of the system 3 . They have 
considered nonrotating beams, and determined the effect of 
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these nonlinearities on the response of the system for the 
cases in which the torsional frequencies of the beam are 
much larger than its bending frequencies. For such cases, 
the nonlinearities present in the differential equations of 
motion are 0(e 3 ) rather than 0(e 2 ). 

The question that immediately arises for the rotating 
rotor-blade problem is whether cubic nonlinearities can also 
play a significant role in the equilibrium and stability of the 
elastic motion of the blade. To address this question, the 
differential equations, and their boundary conditions, de- 
scribing the flap-lead/ lag-torsional elastic motion of a rotor 
blade were derived in Ref. 4 with the objective of retaining 
in the equations all the nonlinear terms up to 0(e 3 ) in a 
small parameter e. The equations developed in Ref. 4 ex- 
tend those developed in Ref. 1 to include not only all linear 
0(e 3 ) terms but all nonlinear terms to this same order. 

In this paper, the 0(c 3 ) differential equations devel- 
oped in Ref. 4 are used to investigate the influence of these 
higher-order terms in the elastic response and stability of 
a rotor blade in the hover flight condition. First, a brief 
review of the derivation of the equations is given. A small 
arbitrary ordering-parameter c is then introduced and the 
equations are simplified by expanding their nonlinearities 
into a power series in e. The resulting equations are more 
amenable to analysis, and Galerkin’s method is applied to 
them. After the equilibrium solution to the equations is de- 
termined, the blade’s elastic deflections are then perturbed 
about their equilibrium to yield a set of variational equa- 
tions that are linearized and used to determine the eigen- 
values associated with the perturbed motion. The influence 
of a number of 0(e 3 ) terms on the blade’s response is de- 
termined for a range of values of collective pitch. 


2. Equations of Motion 

2.1 Basic Assumptions and Outline of Derivation 

Consider an initially straight rotor blade of closed cross 
section. Its maximum cross-sectional dimension is assumed 
to be much smaller than its undeformed length R, so that 
it may be approximated as a beam. A blade segment, both 
in its undeformed and deformed states, is shown in Fig. 1. 
The (t 7 ,f, £) axes shown in the figure, with unit vectors 
indicated by a hat as () , are the principal axes of the cross 
section at the shear center C e of the deformed blade cross 
section. It is assumed that the cross section is symmetric 
about the r/- axis. The £-axis is tangent at all times to the 
elastic axis of the blade. When the blade is undeformed, the 
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Fig. 1 Undeformed and deformed blade segment with co- 
ordinate systems and unit vectors. 


ing axes shown in Figs. 1 and 2. Here, (ti,v,u;) are the 
components of the elastic displacement vector for C e , nor- 
malized by R. They are functions of the nondimensional 
variable x-the distance along the x-direction, normalized 
by R - and of the nondimensional time r. The orientation 
of the cross-sectional principal axes £) centered at G e 
may be described by a sequence of three-axes Euler angles 
(8 Z — arctan v'/(l + u'), 8 y — arcsin iw'dx/dr, 8 X ), as de- 
scribed in Refs. 1-4. Here, primes are used to denote partial 
differentiation with respect to x, and 

§^ = [(1 + «') 2 + «/ 2 + «/’]“ 4 ( 2 . 1 ) 

The elastic angle of twist of the blade, <£(x,r), is ob- 
tained by integrating the torsion of the blade and is related 
to the Euler angle & x (x,t) as 1-4 


principal axes make an angle 0 c -the collective pitch 

angle-with the y-axis. The ( x,y,z ) axes, with unit vectors 
also indicated by a hat , are a set of rotating reference axes 
as shown in Fig. 2. The x-axis is coincident with the elastic 
axis of the blade when it is undeformed. These axes are 
assumed to rotate in space with constant angular velocity 
H about the vertical, which is taken to be perpendicular 
to the rotor hub. The (X,Y,Z) axes shown in Fig. 2 are 
a set of inertial axes. The absolute orientation of (x,y,z) 
may be described by first aligning (x,y,z) with (X, F, Z) 
and then performing two successive rotations. The first 
rotation r = fit, where t denotes dimensional time, about Z 
brings the (x, y, z) triad to its new orientation (X 1} Fi, Z \ = 
Z); a second rotation — the blade’s pre-cone angle — about 
the negative Y\ direction brings (X\yY\yZ\) to its “final” 
orientation (x,y, 2 ). For simplicity, the blade-root offset t\ 
shown in Fig. 2 is assumed to be zero. 

Because of the elastic deformations, point C e in 
Fig. 1 moves from location (i?x, y=0, z= 0) to \Rx + 
Ru(x,t),Rv(x,t),Rw(x,t)} relative to the (x,y, z) rotat- 



Fig. 2 Nonrotating and rotating coordinate systems with 
unit vectors. 


<f> = $ x + f 0^sin 8ydx (2-2) 

J o'* 

To obtain the differential equations of motion, and 
their boundary conditions, use is made of Hamilton’s ex- 
tended principle 5 . These equations were developed in Ref. 4 
in terms of the elastic deformations u(x,r), v(x,r), and 
tt>(x, r),and of the angle 0 z (x,r). If the blade’s mass cen- 
troid offset from its elastic axis is neglected, for simplicity, 
the equations associated with the virtual displacements <5u, 
6v } and 6w are of the form 


Gu(x,r) = u - 2 v cos 0 -j- tu(sin2/?)/2 
-(x + u) cos 2 /? — Q u 

G' v (x,r) = v + 2 ticos /3 — 2 ws\n/3 - v - Q v (2.4) 


G[y(x,r) =w + 2 vsm/3 + (x + u)(sin2/?)/2 
— u>sin 2 /? — Q w 


with the cantilever boundary conditions 


u(0,r) = t/(0, t) = w(0,r) 

= 8 x (Q)T) = t/(0,r) = u/(0,r) = 0 


(2.6) 


G u (l,r) = Gy(l,r) = G w { l,r) = *'(l,r) = 0 (2.7) 


In the above equations, dots denote differentiation with 
respect to t. The G u , G y , G w and G$ z terms are nonlin- 
ear functions of the elastic deformations and of their spa- 
tial and temporal derivatives 4 . The Q u , Q v , Q w terms 
are the distributed forces (normalized by mRO 2 , where 
m is the blade’s mass per unit length, which is assumed 
to be constant) associated with the virtual displacements 
RSuy RSvy and R6w, respectively; and Q$ t is the dis- 
tributed moment, normalized by mi? 2 n 2 , associated with 
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the virtual rotation 66 x . The normalized virtual work 
due to these generalized forces is expressed in the form 
Qu + + Q w 6w + Q$ z 66 x . The boundary condi- 

tions <7 w (1,t) - 0 ^( 1 , r) = 0 imply t/'(l,r) = w"(l,r) = 

t/ w (l,r) = itf w (l,r) = 0. 

For compactness, the fourth differential equation ob- 
tained from Hamilton’s principle, namely, the equation as- 
sociated with the virtual rotation 66 x , is presented in the 
next section in its simplified expanded form only. For its 
complete form, the reader is referred to Ref. 4. 

The aerodynamic forces and moments acting on the 
blade are modeled using quasi-steady strip theory based 
on Greenberg’s extension of Theodorsen’s theory in which 
only the (f 7 ,<r) components of the blade’s elastic axis ve- 
locity relative to the air are assumed to affect the aero- 
dynamic loading 2,4 ’ 6-8 . These components, normalized by 
the blade’s tip speed HR, are given as 4 


U n = T 22 [ v + (x + u) cos ft — w sin ft] 
+T23 [ ti; + v sin ft + A cos ft] 

= r 32 [ v + (x + u) cos ft — w sin ft] 
+T33 [ w + v sin ft + A cos ft] 


( 2 . 8 ) 

(2.9) 


F ” ~ 6 


U 1 ~ ~ + U ! 


(2.12 - c) 


F s = - 
f 6 


- Ur,U s + -U n w ( --u f + -& e 

-% U *'J U 3 + U ? 


( 2.12 - d) 


In the above equations, c denotes the blade’s chord, 
normalized by 72, 7 is the Lock number, c^o is the airfoil 
profile drag coefficient, and is the £ component of the ab- 
solute angular velocity of the principal axis system (r?, f, £). 
It is given as 


= 6 X + ( 6 Z + cos ft) sin 0 y + sin/? cos 6 y cos 6 Z (2.13) 

In order to compare results with those obtained via 
the equations developed in Ref. 2 , the normalized induced 
inflow A is modeled as being uniform along the blade radius 
and is given as 


where, 

T 22 — cos 6 X cos 6 Z — sin 6 X sin 6 y sin 6 Z (2.10 — a) 

T 23 = sin0 x cos0 y (2.10-6) 

T 32 = — (sin 6 X cos 6 Z + cos 6 X sin 6 y sin 6 Z ) (2.10 - c) 

T 33 = cos 6 X cos 6 y (2.10 - d) 

and A is the induced inflow velocity normalized by UR. 

As shown in Ref. 4, the generalized forces Q u , Q v , Q wy 
and Qe x due to the aerodynamic loading are determined as 

Qu = ^ 21 ^ + r3 1 F f -Q e X§^ , ( 2 . 11 -a) 


Qv = T^F^ + T 32 F^ — Qe x w 


, d6 z dx 
du' dr 


( 2.11 - 6 ) 


Qw = t 23 f„ + t 33 f< 


( 2 . 11 -c) 


A = sgn [ 6 C + <£c(0.75)] 


be 




1 + ^ I + ^.(Q- 75 ) I 


-1 


(2.14) 


where b is the number of blades, and <^(0.75) is the equi- 
librium value of the elastic angle of twist at x = 0 . 75 . 


2.2 Ordering Scheme and Expansion of the Equa- 
tions to 0 (e 3 ) 

Because of the complexity of the differential equations 
presented in the previous section, they will now be re- 
stricted to moderately large deflections by expanding their 
nonlinear terms in a Taylor series in a small ordering pa- 
rameter e, and truncating the result to 0(e 3 ). Our ob- 
jective here is to evaluate the influence of these higher- 
order terms on the motion of the system. We then let 
v{x,t) — O(e), ty(x,r) = O(e) and 0 x (x,r) = O(e). In 
addition, u(x,r) = 0(e 2 ) . As an example, the expanded 
form of 6 y = arcsin w’dx/dr is 


Qe x = - 


with, 


16 


(“(\Juz+u?-u i ) + 3 ^ u ( 


( 2.11 - d) 


T 21 = — (cos 6 X sin 0 Z +sin 0 I sin(? !/ cos 0 *) ( 2.12 — a) 


e y - u/(l - u' - t/ 2 /2) - ui' 3 /3 + 0(e 5 ) (2.15) 

By making use of the boundary condition G u (l, t) = 0, 
Eq. (2.3) may be integrated over x to obtain an expression 
for u' in terms of the remaining variables. With u(0, r) = 0 , 
the following expression is obtained for u(x,r) (Ref. 4 ) 


2 ai = (sin 6 X sin 6 Z — cos 6 X sin 6 y cos 6 Z ) ( 2.12 - 6 ) 
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The blade’s normalized area moment, 7^, is 


u(x,r) 


\ / W 2 + w ' 2 ~ Ja t 1 _ cos2 ft dx 

-Jr: f f [2 vcos /? + Q u ] dx dx + 0(e 4 ) 
EAJ 0 Ji 


= 0(e 2 ) 


U = J f (I 2 + $ 2 )dn d< (2.18 -d) 

and (Dfj, D^D{) are the blade’s flexural and torsional nor- 
malized stiffnesses determined as 


where E is the Young’s modulus of the material, normalized 
by mH 2 , and A is the blade’s cross sectional area, normal- 
ized by 72 2 . Both of these quantities are assumed here to 
be constant. 

With u as given by Eq. (2.16), the 0(e 3 ) expansions for 
the quantities G v (x,r) and G w (z,r) in the Sv and Sw equa- 
tions, Eqs. (2.4) and (2.5), respectively, are now in integro- 
differential form. Furthermore, since G v (l, r) = C7 w (l,r) = 
0, it is convenient to reduce these equations further by in- 
tegrating them in x from x = 1 to x = x and applying 
the Galerkin procedure to the latter equations. For sim- 
plicity, it is also assumed that sin/3 = O(e) (with cos/? left 
as an O(l) quantity in the equations), and that c = O(e), 
A = 0(e), and Cdo = 0(e 2 ) in the generalized aerodynamic 
forces. 

The expanded 0(e 3 ) form of the fourth differential 
equation obtained from Hamilton’s principle namely, the 
equation associated with the virtual rotation 6B Xi becomes, 
after some higher-order cross-sectional integrals are ne- 
glected, as is commonly done in the literature (e.g., Refs. 1- 

4 ), 


H (f x + cos ft ) + Or ~ Jv){ v f sin(2fl c ) - u>'cos(2fl c ) 

+ [O x cos(20 c ) + sin0 c cos B c \ cos/?} cos/? 

+ ( D n - ( « " 2 - w " 2 ) sin2(fl c + 0 X ) 

- v"w ff cos 2 (B c + B x )\ 

+ -— 7 - [ x ^ 2 0 Z + w' cos /? + sin/? ^ cos /? cos 6 C 
+ v sin 8 C — w cos B c ] A 0(e 4 ) = 0 

(2.17) 

All quantities in the above equation are nondimen- 
sional. The blade’s distributed mass moments of inertia 
are determined in terms of its material density 

P{ru$) as 


mj n =R 2 J j p?dr\ d{ 

(2.18 - a) 

= Rl J j k 

(2.18 - 6) 

3i~ Jv + 

(2.18 - c) 


D V = E J J s 2 d„ d( ( 2.18 — d) 

D,=eJ J ^dr, ds (2.18 - e) 

D( = G f flk + ^ + il- f :) 2 ]*l dt (2.18 - /) 

where G is the normalized shear modulus of the blade’s 
material and is the warp function (normalized by 

72 2 ) for the blade’s cross section. It is assumed that t/> is 
anti-symmetric in (*?,£)* 

For compactness in presentation, two terms in 
sin2(0 c + 0 z ) and cos2(0 c + 0 x ) are shown in Eq. (2.17), 
but they were actually approximated by their respective ex- 
pansions to 0(e 3 ) about $ x = 0 by writing sin2(0 c + B x ) = 
sin 2 B c + 2 B x cos (2 B c ) + 0(e 2 ) and cos 2{B C + B x ) = cos 2 B c — 
20 x sin(20 c ) + 0(c 2 ). The nonlinear 0(c 3 ) terms associated 
with these expansions, namely, the terms in v ,,2 B x , w ,,2 B x 
and v”w ,, B Xi will henceforth be referred to as the ijkl terms 
in the 6B X equation in the next figures. The 0(e 3 ) terms 
underlined in Eq. (2.17) are not included in the equations 
developed in Ref. 2. The single underlined terms are lin- 
ear pitch-flap and pitch- lead/ lag coupling terms, and the 
remaining underlined terms are 0(e 3 ) linear terms in the 
aerodynamic pitch moment that are kept for consistency 
in the formulation. Until a better understanding and more 
accurate modeling of aerodynamic phenomena is achieved, 
the validity of terms such as these may be questionable. 
The l/(EA), 0(e 3 ), terms in Eq. (2.17) were also neglected 
in Ref. 2. Again, these terms are kept here for mathemati- 
cal consistency. For values of EA greater than about 200, 
we found that the influence of these terms in the results pre- 
sented later is so small that they may actually be neglected 
in practice. 

2.3 Application of Galerkin’s Method 

We approximate the solution to Eq. (2.17) and to the 
integrated form of Eqs. (2.4) and (2.5) as a series of the 
form 

N 

v(x,r) = ^2 v tj(r) fj{x) (2.19-o) 

j=i 

N 

w(x,r) = 52 w ‘j( t ) fj( x ) (2.19 - b) 

> = i 

N 

M*,r) = H M r ) 9ii x ) (2.19 - c) 

>=» 
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and then reduce the integro-differential equations to or- 
dinary differential equations by making use of Galerkin ’s 
method Ref. 5. The functions /y(x) and gj(x) are cho- 
sen here as the orthogonal eigenfunctions for a nonrotating 
clamped/free beam, 


where 


Ri 


i / X y,' d,2 

Jo 

= ( l x 9i fj 
Jo 


dx 


(2.22 - a) 

(2.22 - 6) 


/y(x) - cosh (0jx) - cos (Pjx) 

- ay [sinh (Pjx) - sin (Pjx) J 

(2.20 - a) 

L 3,i} = - fi'gj dx 

Jo 

(2.22 - c) 

9 j(x) = N/2sin[7r(j - §)x] 

(2.20 - 6) 

°4,ij = f fi g, dx 

Jo 

(2.22 - d) 

where 


B3,ijki = / fi'gjgkfi" dx 

Jo 

(2.22 - e) 

cos 0.- + cosh 

“’~.i»0,+.m hft I 2 - 20 -') 

and is the (j = 1,2 N) root of the characteristic 

equation 

The terms that are underlined in Eq. (2.21) correspond 
to those similarly underlined in Eq. (2.17). The L 3| ,y, 0 4>1 y, 
and Bz,ijki Galerkin coefficients also appear in the 6v and 
in the 6w equations, with the I 3|i y coefficient in the form 

of a $tj term. The S+ and i2,y coefficients appear only in 
Eq. (2.21). 

1 + ( cosh 0j ) cos Pj = 0 

(2.20 - d) 




All the Galerkin coefficients obtained by the proce- 
dure described above were evaluated numerically, stored 
in a computer file, and then used to generate the results 
presented in Sections 3 and 4. The ordinary differential 
equation obtained by applying Galerkin’s method to the 
86 x equation is obtained as 


2 ^ -^)Qi( s in2^ c )cos 2 ^ + ^(sin2/?)(cos 0 C )J^ 

N / ~ 

+ ( [■***«> + C?f “ cos 2 0 C ) cos 2 0]6ij 

3 - 1 v 

” + Of — Jv){vtj sin 2 9 C — w t j cos 2 0 C )] L 3> y, cos 0 


+ ( cos P) (cos 0 C ) [2 0 t; M,y 

+ (cos £)w«yjfry] + ( V tj sin 0 C - tt> t y COS g c ) P 4| y, | j 
JV JV , 

+ EE DiMjikVtjWtk 

3 = i *=i 1 

+ (- 0 *? - ^f)[ ~ wtjtvtk) sm20 c 


VijWtk cos2$ c \K, 


i;fc| 


N N N 

+ (A? ~ E)f ) ^ (Vt; Vtfc - COS 20 c 

;'=i *=i /-i 

+ 2 sin 20 c ]0 t/ J9 3> y l/jt 

+ 0(e 4 ) = 0 (t = 1,2, ...,7V) 


( 2 . 21 ) 


3. Equilibrium Solution 
3.1 Numerical Method 

The differential equations outlined in Section 2 admit 
the equilibrium solution v t y(r) = constant = v e y, w tj {T) = 
constant = u; e y, and 0 f y(r) = constant = 0 e y. The 37V quan- 
tities v e y, t v ej and 0 c y, j = 1,2, ...,7V, were determined nu- 
merically by solving the algebraic equations obtained from 
the differential equations in Section 2 using a minimization 
program 9 . 

The equilibrium solutions were obtained for a four- 
bladed rotor with c = 7r/40, and using a Lock number 7 = 5, 
a profile drag coefficient c d0 = 0.01, and EA = 200. The 
equilibrium deformations at the blade tip, v c (x = l),tu c (x = 
1) and (f> e (x = l) are plotted in Figs. 3 and 4 versus 6 C for 
P ~ 0. The quantities u>*, wjj,, and cjJ shown in these and in 
the subsequent figures denote, respectively, the first rotat- 
ing uncoupled blade natural frequencies normalized by 0 
as obtained in Ref. 10. The results shown in these figures 
were obtained by using TV = 5 nonrotating beam normal 
modes in the Galerkin procedure. Greater values of TV did 
not significantly affect the results obtained. The dashed 
lines shown in Figs. 3 and 4 and in subsequent figures rep- 
resent the results obtained using the equations in Ref. 2, 
while the solid lines represent the results obtained when 
the additional 0(e 3 ) terms presented here are included in 
the differential equations of motion. These lines are marked 
a, b, c and d and they represent the following cases: 

a the full 0(e 3 ) equations; 

b the 0(e 3 ) equations, but with B Z} ijki = 0; 

c the 0(c 3 ) equations, but with all the 0 4(t y, R t; and 
Lz y ij terms removed, and Bz,ijki = 0; 

d the 0(c 3 ) equations, with all the 0 4> ,y and R i: terms 
removed. 
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3.2 Discussion of Results 

Figs. 3 and 4 illustrate the importance of the 0(c 3 ) 
terms of the type underlined in Eq. (2.21) in the equilib- 
rium solution of the blade. Their effect is particularly re- 
flected in the equilibrium value of the blade’s angle of twist, 
<f> e (x = 1), especially for the lower value of the uncoupled 
flap frequency = 1.06. Several of the nonlinear terms 
may have only a minor influence on the blade’s equilib- 
rium deflections. For soft in-plane blades with = 0.7, 
for example, it is seen that the equilibrium curves a and 
b are nearly identical; the same is true of curves c and d. 
This indicates that for such blades, the additional aerody- 
namic 0(e 3 ) terms of the type indicated in Eq. (2.21) are 
essentially responsible for the difference between curve a 
and the remaining curves. For these blades, the nonlin- 
ear B 3t ijki terms in Eq. (2.21) could have been neglected 
without causing any significant change in the blade’s equi- 
librium deflections. For high values of collective pitch, how- 
ever, those terms exert a significant influence on the blade’s 
equilibrium, affecting especially its elastic angle of twist. 

Another characteristic of the full 0(e 3 ) equations is 
disclosed by examining the equilibrium response of a stiff 
in-plane blade with wj, = 1.06. The numerical determina- 
tion of the equilibrium deflections based on the equations 
developed in Ref. 2 fails to converge when 9 C is about 0.4. 
This singularity is shifted to a higher value of 9 C when ad- 
ditional 0(e 3 ) terms are included in the equations. For the 
full 0(e 3 ) nonlinear equations, no singularity is exhibited in 
the range of 0 C shown in Figs. 3 and 4. If the aerodynamic 
04,»y and R{j terms and the Bz^jki terms are neglected, 
but if all other additional terms in the equations are kept, 
the singularity now appears near 0 C = 0.5. 

It was verified that the 0(t 3 ) ijkl terms in the 6v and 
Sw equations , did not contribute significantly to the deter- 
mination of the blade’s equilibrium response and, therefore, 
could have been neglected for practical purposes. 


4. Stability Analysis 
4.1 Numerical Method 

To analyze the stability of the motion about the equi- 
librium determined in Section 3, we let 

v tji T ) = t>ey+t>/(r) (4.1 -a) 

w t j{r) = w e j + tuy(r) (4.1 - 6) 

*t/( r ) = Qej + 0/(r) (4.1 - c) 

and then linearize the 3 N differential equations of motion 
in the variables vy(r), UJ ; (r), and 0y(r) to obtain a matrix 
equation of the form 

M‘q + C £ + Kq = 0 (4.2) 

where £ is a 3 N x 1 column vector whose components are 
vj(r), ... , VAr(r), ^(r), ... , u; N (r), ^(r), ... , 6 n (t). 
The matrix M is symmetric, and the matrices K and C 
are non-symmetric. 


The stability of the perturbed motion g(r) is deter- 
mined by the eigenvalues associated with Eq. (4.2). To 
determine such eigenvalues, Eq. (4.2) is first rewritten in a 
first-order form. After introduction of a column vector z 
with components q and £, Eq. (4.2) may be written as 

B t = Az (4.3) 

with Bn = 7, a 3 N x 3 N identity matrix; B\ 2 — B 2 \ = [0]; 
a 3N x SN null matrix; B 2 2 = M; An = [0]; A 12 = /; 
^21 = and A 22 = —C. The eigenvalues associated 
with the 6 N x 6iV matrix in Eq. (4.3) were determined 
numerically by making use of the IMSL routine EIGZF®. 

The real and imaginary parts of the first lead/lag ( a v 
and Wt,), first flap (o w and w w ), and first torsion ( o $ and 
u<f>) eigenvalues determined as indicated above are plotted 
versus collective pitch (0 C ) in Figs. 5 to 10 using the same 
parameter values and labeling convention indicated in Sec- 
tion 3.1. 


4.2 Discussion of Results 

Figs. 5 and 6 show the first lead/ lag eigenvalue asso- 
ciated with Eq. (4.3) for the rotating blade as a function 
of the pitch angle B c . It is seen that for a soft in-plane 
blade with uncoupled rotating natural frequency u;* = 0.7 
the nonlinear 0(e 3 ) B 3fi j k i term that appear in Eq. (2.21) 
has no substantial influence on either a v or u v . For such 
blades, the influence of the remaining higher-order terms 
underlined in Eq. (2.21) is reflected in the difference be- 
tween curve d (i.e., with the aerodynamic terms 0 4 t y and 
R\j removed from the equations) and that obtained with 
the full 0(e 3 ) equations (curve a), and curves a and c (i.e., 
with the £3,1; torsion-bending coupling terms, and the 04,»y 
and R{j aerodynamic terms removed from the full 0(c 3 ) 
equations). The additional 0(e 3 ) terms included in the 6v 
and equations account for the difference between the re- 
sults obtained by using the equations developed in Ref. 2 — 
represented by the dashed line — and those represented by 
curves a in Figs. 5 and 6. 

As uj* is increased, however, the situation described 
above changes. For a stiff in-plane blade with u* = 1.5, the 
nonlinear B 3 ,ijki term that appears in Eq. (2.21) now exerts 
a major influence on the real part a v of the first rotating 
lead/lag eigenvalue, whereas the terms in the underlined 
coefficients in Eq. (2.21) do not. The effect of the 0(e 3 ) 
nonlinearities that appear in the 6v and Sw equations is 
seen by comparing curve c with the dashed curve obtained 
by using the equations in Ref. 2. For values of 0 C as high 
as about 0.4, the latter equations yield, for this stiff in- 
plane blade, practically the same values for as the full 
0(c 3 ) equations used in this paper. At about 0 C = 0.4, 
the numerical calculation of the eigenvalues based on the 
equations in Ref. 2 fail to converge. 

The first flap eigenvalue obtained from Eq. (4.3) is 
shown in Figs. 7 and 8. Again, the effect of B 3tijkl on 
both o w and w w is negligible for a soft in-plane blade with 
u;* = 0.7, but significant for a stiff in-plane blade with 
w* = 1.5 for higher values of collective pitch. The effect 
of the apparent inertia, 04 jt y , and of the R^j aerodynamic 
terms, generally neglected in the literature, is reflected in 
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the eigenvalue (ffu,,^) even for 0 C = 0- It was verified 
numerically that the Lz t ij torsion-bending coupling term 
that appears in Eq. (2.21) — and also in the Sv and 6w vari- 
ational equations as a term — has no practical influence 

on tJ w . Its influence on o w was found to be negligible for 
u* = i.5 ? but significant when u* = 0.7. As indicated by 
Figs. 7 and 8, the damping for the perturbed flap motion 
for a soft in-plane blade can be significantly affected by the 
additional 0(c 3 ) terms included in the equations used here. 

co* v = 0.7 




Figs. 9 and 10 show the real and imaginary parts of the 
first rotating torsion eigenvalues of Eq. (4.3). The values 
of the torsional frequency are relatively large, and, as 
seen from these figures, there is little difference between 
the results obtained here and those obtained by using the 
equations in Ref. 2 for small values of 0 C . For larger values 
of collective pitch 0 C , however, the full 0(c 3 ) equations used 
here and those in Ref. 2 predict a different trend for as 
6 C is increased further. However, this trend difference is 




Fig. 9 First torsion eigenvalues (real part, a^, and imaginary part, oo 0 > versus collective pitch (0 C ) for 
oo* = l.O6,0 = O,o;; = 2.5,(N = 5). 





Fig. 10 First torsion eigenvalues (real part, o^, and imaginary part, to^) versus collective pitch (0 C ) for 
oo^ = 1.15,0-0, = 2.5, (N = 5). 


exhibited at such high values of B c that it may not be of 
practical significance. 

The torsion damping is significantly affected by the 
0(e 3 ) terms included in the equations used in this paper. 
Particularly noticeable in Figs. 9 and 10 is the opposing tor- 
sion damping trends for a soft in-plane blade with w* =0.7 
for increasing values of 6 C . This is observed by comparing 
the results given by curve a, obtained with the full 0(e 3 ) 
equations, and by the dashed curve obtained by using the 
equations in Ref. 2. 

It is worth mentioning that all the results shown in 
Figs. 3 to 10 were essentially unaffected by the 0(e 3 ) ijkl 
terms in the Sv and Sw equations . These are terms in 
vtjtvtkOth VtjWkOth etc. 


5. Concluding Remarks 

Numerical results obtained from nonlinear rotor blade 
equations for the hovering flight condition, with terms re- 
tained up through 0(e 3 ), are presented and compared with 
results from a simpler 0(e 2 ) model obtained by previous 
investigators 2 . In order to facilitate an understanding of 
which terms are important in the present model that were 
absent in the previous, simpler model, the present model 
was exercised with several different classes of terms sys- 
tematically omitted. Present results, a subset of all the 
results obtained, indicate that both linear and nonlinear 
terms of 0(e 3 ) can significantly affect results for both non- 
linear static equilibrium and linear aeroelastic stability. For 
the results presented here, the most significant cubic non- 
linear terms are those associated with structural geometric 
nonlinearity in the torsion equation. It would appear that 
such terms should be present in any general-purpose rotor 
dynamics analysis. The corresponding terms in the equa- 
tions for bending, although not practically significant in the 
present results, do make the structural terms in the equa- 
tions symmetric. The most significant linear terms in the 
present model but absent in Ref. 2 are associated with an 


approximate aerodynamics model, the accuracy of which 
has not been rigorously ascertained. For completeness, it is 
recommended that a similar investigation be undertaken for 
the forward flight condition to determine if similar trends 
hold. 
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